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ON THE COMPUTATION OF SOLAR ELEVATION ANGLES 
AND THE DETERMINATION OF SUNRISE AND SUNSET TIMES 


by 

Harold M. Woolf 


National Meteorological Center, ESSA 
Hillcrest Heights, Maryland 

SUMMARY 


The complete procedure for precise computation of solar elevation angle as 
a function of latitude, longitude, date, and time is given. Construction of a 
graphical aid for determining times of sunrise and sunset, with a precision of 
one minute, as functions of latitude, longitude, date, and altitude is described. 
This Sunrise- Sunset Finder is contained in the Appendix. 
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ON THE COMPUTATION OF SOLAR ELEVATION ANGLES 


AND THE DETERMINATION OF SUNRISE AND SUNSET TIMES 


by 

Harold M. Woolf 

Environmental Science Services Administration 
Weather Bureau 

National Meteorological Center 


INTRODUCTION 

Upper -atmosphere research workers often require solar elevation (or zenith) 
angles as input for various problems. Examples include the determination of 
solar -radiation corrections* for stratospheric radiosonde data and the calcu- 
lation of photo-dissociation and ionization rates in the mesosphere and ther- 
mosphere [ 2 ]. Solar elevation as a function of declination, local time, and 
latitude is presented in convenient graphical form in the Smithsonian Meteor- 
ological Tables [3]. However, for many purposes the charts do not permit solar 
elevation to be determined with sufficient precision. In addition, many of the 
calculations in which these angles are employed are performed by electronic 
computer. Thus it seems worthwhile to summarize the steps required for precise 
determination of solar elevation. The procedure is quite readily programmed 
for computer operation; one such program has been in daily use sinee the 
beginning of 1964 [4]. 

On the other hand, there are applications of earth-sun relationships in 
which direct calculation is not the optimum approach. For example, it is often 
necessary in connection with the planning of such experiments as chemical trail 
releases from rockets, to know the time of sunrise and/or sunset at a particular 
location and altitude. While tables of sunrise/sunset times at the earth's 
surface are readily available [ 5 ], adjustment for altitude entails tedious 
exercises in earth- sun geometry and computations with trigonometric functions 
as in the case of solar elevation. However, with proper manipulation of the 
relevant astronomical parameters and formulas, it is possible to perform the 

^These corrections are based on the average differences between daytime and 
nighttime observations, which have been found [lj to depend strongly on the 
solar elevation angle of the daylight sounding. They are needed to achieve 
compatibility of data, obtained at various local times and with many different 
instruments, for purposes of synoptic analysis [4], 
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Figure 1. Solar declination (D) and time of meridian passage (M) versus day of the year. 


tedious calculations once and for all, and construct tables or graphs to 
expedite and simplify the determination of sunrise/ sunset times for any 
altitude and location. 

CALCULATION OF SOLAR ELEVATION ANGLE 

The solar elevation angle o( is determined by the relation 

sin o( = sine/ sin D + cos <j> cos D cos h (.1.1) 

where <fi is the latitude, D the solar declination, and h the solar hour angle. 
At first glance, declination (Fig. l) is a simple sinusoidal function of date, 
with maximum and minimum at the summer and winter solstices, respectively, and 
nodes at the equinoxes. Careful inspection of the curve, however, reveals a 
slight asymmetry, due to the ellipticity of the earth's orbit, which is 
accounted for in the following procedure for computing declination*: 

sin D = (sin 23°26'37.8") sin O' (1.2) 

where cr( deg) = /+ 0.4o87sin./ + 1.8724cos>/ 

- 0.0182 sin 2.1 + 0 . 0083 COS zl (l- 3 a) 

i(deg) = 279.93^8 + d, (1.3b) 

and d is the angular fraction of a year represented by a particular date, 
given by 

[(number of day in year)-l]x360 deg. For example, on 1 January, d=0 
355.21+2 

and on 21 March, d = 79(0.98565) - 77-86635° = 77°52'. 

It is convenient to combine equations (l.3a) and (l. 3b) to yield 

CT = 279.93W + d + 1.9l4827sin d - O.O79525COS d 

+ 0.019938sin 2d - 0.00l620cos 2d (l.U) 

The solar hour angle h, a measure of the longitudinal distance to the sun 
from the point for which the calculation is being made, is given by 

h(deg) = 15(T-M)-L (l-5) 

♦Equations (12) and (l.3) and the definition of d used here were obtained from 
the Nautical Almanac Office, U.S. Naval Observatory, Washington, D.C. Year- 
to-year variations in the relationships involving <r, £, , and d are negligible 
for atmospheric applications. 



where T(hr) is the time (GMT) of the calculation! M(hr) is the time of meridian 
passage, or true solar noonj and L(deg) is longitude, counted positive west of 
Greenwich. M, which is shown along with D in Fig. 1, is given by 

M(hr) = 12 + 0.123570sin d - 0.004289cos d 

+ 0.153809sin 2d + 0.060783cos 2d (1.6) 

where d is as defined earlier. 

In the operational program for applying radiation corrections to high- 
altitude radiosonde data [4], T in equation (1.5) is replaced by H Q + c, 
where H Q is a standard synoptic observation time such as 0000 or 1200 GMT, and 
c accounts for the time elapsed between release of the balloon and its esti- 
mated time of arrival at successive standard pressure levels. Since actual 
release time and ascent rate are not included in upper -air reports transmitted 
over the meteorological telecommunications network, representative values of 
these parameters along with standard- atmosphere [6] heights have been employed 
to construct a table of c-values for use in the program (Table l). 


TABLE 1. Corrections to nominal observation time 
based on estimated time of radiosonde arrival at 
selected constant-pressure surfaces using ascent 
rate of 305 meter s/minute (after [4]). 


P(rnb) Z(m) 


c e (hr) 


C£(hr) 


100 

l 6 l 80 

70 

18442 

50 

20576 

30 

23849 

20 

26481 

10 

31055 

7 

33453 

5 

35791 


0.134 

0.551 

0.258 

0.675 

0.374 

0.791 

0.553 

O .970 

O .697 

1.114 

0.947 

1.364 

1.078 

1.495 

1.206 

1.623 


c g : for "early release," nominally 45 min. prior to observation time. 

C£! for "late release," nominally 20 min. prior to observation time. 



A GRAPHICAL AID FOR DETERMINING SUNRISE AND SUNSET TIMES 


Derivation 


The solar day extends from M-H to M+H, where M is the time of meridian 
passage and H is the solar hour angle corresponding to the time of the sun’s 
appearance or disappearance on the horizon. The latter phenomena are described 
uniquely by the sunrise/ sunset solar elevation angle A, which is measured with 
respect to the earth's surface and therefore decreases (increases in the 
negative sense) with height. Once the variation of A with altitude is known, 

H can be found for any height, latitude, and solar declination. 

From the earth- sun geometry, values of A were determined at 5 -km intervals 
from the surface to 120 km, and corrected for refraction with an assumed solar 
angular radius (limb to center of disk) of 0.25 deg. The variation of A(deg) 
with geometric height Z (km), shown in Fig. 2, is described almost exactly 
(mean error -0.005 deg, r.m. s. error 0.057 deg) by the relation 

A = -1.76459Z°' 1K)795 (Z>l), (2.1) 

obtained by least squares. 

Equation (l.l) may be rewritten, for the special case of sunrise/sunset, as 

sin A = sin^sin D + cos</ cos D cos H (2.2a) 

„ f sin A- sin <±> sin D / 

or H = arecos — — -r-*- > x 

I cos <p cos D ( (2.2b) 

Equation (2.2b) was solved by high-speed computer with A varying from -0.9° 
(surface) to -12.5° (120 km) in 0.1-deg increments, from 0 to 85° in 1-deg 
increments, and D from - 24° to 24° in 1-deg increments. The resulting values 
of H were printed out in units of hours and minutes. 


Presentation 


After considerable experimentation, the following method was chosen as 
optimum for presenting the results of the above computations in a convenient 
and readily usable form. For each of 10 latitudes (0, 10, 20, 30, 40, 4-5, 50, 
55> 60, and 65°) a chart was constructed on which H is a function of height for 
a given declination. The resolution in declination on individual charts, and 
in latitude between charts, is such that values of H obtained by linear 
interpolation are accurate to one minute. Graphical display is not practical 
for latitudes above 65°, where H changes very rapidly with declination, 
especially near the summer solstice. 

The complete Sunrise-Sunset Finder is contained in the Appendix. In order 
to avoid excessive photo-reduction and consequent loss of readability, the 
charts for 55> 60, and 65° are each presented in two sections. 




A(DEG) 


Figure 2. Variation with height of solar elevation angle, relative to the 
earth’s surface, at sunrise and sunset. Equation was obtained by least- 
squares fit to points at 5 -km intervals. 
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Once H has been found, times of sunrise and sunset are given in Greenwich 
Mean Time (GMT) by the relations 

t sr - 15 + M - H ( a - 3s > 

T SS " T5 + M + H (2.3b) 

a where L is longitude in degrees, counted positive west of Greenwich. All 

quantities in equations (2.3a, b) are in units of hours and minutes. If the 

time zone of the point is known, Tg^ and Tgg may readily be converted to 
Local Standard Time if desired. 
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APPENDIX 


SUNRISE- SUNSET FINDER 
INSTRUCTIONS 

* 

The Finder permits the determination of times of sunrise and sunset for any 
altitude from sea level to 120 km, any latitude up to 65°, and any longitude 
and date . a 

Example : 

Determine times of sunrise and sunset at 55 km, at 39°R 77°W, on 
January 1. 

(1) From Fig. 1 of the main text, extract solar declination D (degrees) 
and time of meridian passage M( hours and minutes) for the desired date. 

D - -23°, M = 12 hr 03 min. 

(2) Charts in this Appendix are plots of hour angle H at sunrise and 
sunset versus height at selected latitudes, with declination as the defining 
parameter for each curve. Select the chart for the latitude closest to that 
desired, and extract the hour angle for the desired height using linear 
interpolation in declination. Repeat for the chart closest in latitude on the 
other side of that desired, and obtain the required hour angle by linear 
interpolation in latitude. 

latitude 40°; H = 05 hr 29 min. 

30°: H = 05 hr 1*9 min. 

39°: H = 05 hr 31 min. 

(3) Divide the longitude L (degrees, positive west of Greenwich) by 15. 

The integral part of the quotient is then in hours; multiply the remainder by 
60 to convert it to minutes. 

15 = S = 5 + T5 = hr 08 min. 

(1*) Determine the times (GMT) of sunrise and sunset in the following 
manner (all quantities are in hours and minutes): 

(a) Sunrise: t sr = f 5 + M - H = 0508 + 1203 - 0531 = 11^0 GMT 

(b) Sunset: T gg = ^ + M + H = 0508 + 1203 + 0531 = 22l*2 GMT 
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